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ABSTRACT 

Using Spitzer Space Telescope IRAC (3.6-8/im) and MIPS (24/im) imaging, as well as Hubble Space Tele- 
scope optical observations, we identify the IRAC counterparts of the luminous power sources residing within 
the two largest and brightest Lyman-a emitting nebulae (LABs) in the SA22 protocluster at z = 3.09 (LAB 1 
and LAB 2). These sources are also both submillimeter galaxies (SMGs). From their rest- frame optical/near- 
infrared colors, we conclude that the SMG in LAB 1 is likely to be starburst dominated and heavily obscured 
(A v ~ 3). In contrast, LAB 2 has excess rest- frame ~2/im emission (over that expected from starlight) and hosts 
a hard X-ray source at the proposed location of the SMG, consistent with the presence of an active galactic nu- 
cleus (AGN). We conclude that LAB 1 and LAB 2 appear to have very different energy sources despite having 
similar Lyman-a spatial extents and luminosities, although it remains unclear whether on-going star-formation 
or periodic AGN heating is responsible for the extended Lyman-a emission. We find that the mid-infrared 
properties of the SMGs lying in LAB 1 and LAB 2 are similar to those of the wider SMG population, and so it 
is possible that extended Lyman-a haloes are a common feature of SMGs in general. 
Subject headings: galaxies: active — galaxies: high-redshift — infrared: galaxies 
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1. INTRODUCTION 

Narrow-band surveys have revealed an intriguing new pop- 
ulation: giant, radio-quiet Lyman-a emission-line nebulae 
(Steidel et al. 2000), which have been termed 'Lyman-a 
Blobs' (LABs). These objects typically have linear extents 
of 10-100 kpc, and Lyman-a luminosities of 10 43-44 ergs -1 , 
with the highest concentration currently known in the S A 22 
proto-cluster at z = 3.09, with 35 confirmed objects (Matsuda 
et al. 2004). Similar objects have been discovered in other re- 
gions associated with rich, but primitive environments (Keel 
et al. 1999; Francis et al. 2001; Palunas et al. 2004), and ex- 
tended Lyman-a emission has been detected around certain 
high-redshift radio galaxies (e.g. Reuland et al. 2003); how- 
ever, the narrow-band selected LABs are typically not associ- 
ated with such radio-loud systems. 

Cooling-flows or shock heating by outflows originating 
from starburst winds or AGN jets have been proposed as 
the most plausible mechanisms for producing these extended 
Lyman-a haloes (Taniguchi & Shioya 2000; Ohyama et al. 
2003), although other processes such as sub-sonic pdV heat- 
ing, and inverse Compton scattering of Sunyaev-Zelodovich 
photons could also play a role (Scharf et al. 2003). The 
two largest known LABs (LAB 1 and LAB 2, with spatial ex- 
tents of 1 . 1 x 10 4 kpc 2 and 0.8 x 10 4 kpc 2 respectively, both at 
z = 3.09 in SA22; Matsuda et al. 2004) are both associated 
with submillimeter galaxies (SMGs) (Chapman et al. 2001). 
Several other LABs in the SA22 structure at z = 3.09 con- 
tain SMGs (Geach et al. 2005) - the apparent association of 
these bolometrically luminous galaxies with LABs provides 
support for the scenario where feedback mechanisms are pow- 
ering the Lyman-a emission, since the bolometric luminosi- 
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ties of the SMGs are several orders of magnitude larger than 
the Lyman-a luminosities of the haloes. The action of super- 
wind feedback has recently been given credence by integral 
field data which reveals outflows and complex velocity struc- 
tures in both LAB 1 and LAB 2 (Bower et al. 2004; Wilman et 
al. 2005). However, the precise location and identification of 
the SMG counterparts, and hence accurate constraints on the 
dominant power sources, are unclear due to the poor spatial 
resolution of the submillimeter data (Chapman et al. 2004). 

In this Letter, we present Spitzer Space Telescope (SST) 
mid-infrared and Hubble Space Telescope (HST) optical 
imaging of the galaxies embedded within LAB 1 and LAB 2 
to locate the bolometrically luminous sources within them - 
the SMGs. To identify these galaxies, we exploit the fact that 
mid-infrared observations have been shown to be effective 
at identifying SMGs, particularly at 8/im since these galax- 
ies have red infrared colors (Ashby et al. 2006). The mid- 
infrared photometry also constrains their rest-frame optical- 
radio spectral energy distributions (SEDs) and hence can shed 
light on what is powering these sources. Throughout we adopt 
Q m = 0.3, Q\ = 0.7 and Ho = 75 km s -1 Mpc -1 . In this geome- 
try, \" corresponds to 7.1 kpc at z = 3.09. 

2. OBSERVATIONS & REDUCTION 

The SST observations used here are part of GTO project 
#64, which we retrieved from the Spitzer Science Center 
(SSC) archive. LAB 1 and LAB 2 are covered by all four 
IRAC channels (3.6-8/im), and we use the Post Basic Cali- 
brated IRAC frames generated by the pipelines at the SSC. 
Both LABs were also observed by MIPS at 24/im as part 
of the same project, and for these data we perform post- 
processing on Basic Calibrated Data frames to remove com- 
mon MIPS artifacts and flatten small and large-scale gradients 
using 'master- flats' generated from the data. For mosaicking 
we use the SSC mopex package, which makes use of the sup- 
plementary calibration files which are supplied with the main 
science set. We use SExtractor (v2.4.4, Bertin & Arnouts 
1996) to detect and extract sources in the mid-infrared im- 
ages. Fluxes are measured in apertures with diameters chosen 
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Fig. 1 .— (top) Greyscale IRAC 40" x 40" images of LAB 1 (left) and LAB 2 (n'g/tf). The main panels show the IRAC 3.6/xm image, and the 
insets show the 4.5/xm 5.8/xm and 8/xm images, all which have been slightly smoothed with a Gaussian kernel. We overlay the narrow-band 
Lyman-a images from Matsuda et al. (2004) to show the extent of each halo. IRAC detections within the Lyman-a haloes are labelled a, b, 
etc. In the field of LAB 1 we indicate the location of the known Lyman-break galaxies (LBGs; C 11 & C 15) and the position of the 21 cm 
radio source in the halo (white cross, Chapman et al. 2004). LAB 2 contains one known LBG (M 14) and an SMG, as well as an X-ray source 
(white cross). In both cases we show the 15" SCUBA beam (dashed circle) to illustrate the potential uncertainty in the position of the submm 
source, (bottom) High resolution HST-STIS and ACS imaging for LAB 1 and LAB 2 showing the optical counterparts, with contours showing 
the 3.6/xm emission. The white 1" circle in the STIS inset shows the location of the ERO detected in Steidel et al. (2000). In all images, north 
is up, east is left and the horizontal bar is a 5" scale. We discuss the identifications further in §3 and Table 1. 



to be the same as in the SWIRE ELAIS-N1 catalog 6 : 3.8" for 
IRAC and 12" for MIPS. 

In addition to the mid-infrared data, we analyse optical 
HST observations of the LABs, to reveal the rest- frame ultra- 
violet (UV) morphologies of the galaxies in the LAB haloes. 
The STIS coverage of LAB 1 has been previously discussed 
by Chapman et al. (2004), while the Advanced Camera for 
Surveys (ACS) F814W (/-band) imaging of LAB 2 is pre- 
sented here for the first time. The ACS data was reduced 
using MULTIDRIZZLE, and will be described in more detail 
in J. E. Geach et al. (2006, in preparation). 

3. IDENTIFICATION OF THE SMGS 

As the precise location of the SMGs within the LAB haloes 
is uncertain, we begin our discussion by identifying the mid- 
infrared counterparts to the SMGs (which we assume to be 
the dominant power-source of the LABs). Starting with a 
catalogue of the IRAC-detections within each Lyman-a halo 7 
(Table 1), we aim to eliminate sources not at z = 3.1, and then 
identify the most likely mid-infrared counterpart (if any) to 

6 see Surace, J. A. et al. 2004, VizieR Online Data Catalog, 2255 

7 For mid-infrared detections within the Lyman- a haloes, we adopt the 
alphabetical suffix LAB la, LAB lb, etc., although this does not imply that 
these objects are physically associated with the LAB. 



the SMG in each case. Note that no sources are detected at 
24/im within the extent of either LAB, and so we adopt 3-a 
upper limits of /24 < 330 /iJy. 

Figure 1 shows HST and IRAC images of the LABs cover- 
ing 3.6-8/im, and we label each of our potential counterparts. 
Of the five IRAC components in LAB 1 (a-e), we can elim- 
inate several of them as the SMG counterpart immediately: 
LAB lc and d have mid-infrared colors /3.6/A5 > 1, making 
them inconsistent with z = 3.1 galaxies (Ashby et al. 2006), 
while LAB le is at the position of the Lyman-break galaxy 
(LBG) C 1 1, and is only weakly detected at 3.6/im and unde- 
tected at 4.5-8/im. Although this LBG is close to LAB 1, it 
appears to be kinematically distinct from the bulk of the halo 
(Bower et al. 2004) and so we discount it - leaving LAB la 
and b as potential counterparts. 

The positional uncertainty between LAB la and the 850/im 
centroid (A# = 2.0 ± 1.9") is consistent with it being the 
source of the submm emission. The remaining IRAC coun- 
terpart to the south-east of the SMG (LAB lb) is just within 
the Lyman-a halo, and has a bright counterpart in the STIS 
imaging. It is not clear from the mid-infrared colors whether 
this galaxy is at the halo redshift, but given its large separation 
from the nominal submillimeter position, it is unlikely that it 
is associated with the power-source. However, if LAB lb is at 
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Fig. 2. — Rest-frame ultraviolet to radio spectral energy distributions (SEDs) of the proposed power sources in LAB 1 and LAB 2. We overlay 
the SEDs of M 82 and Arp 220 (Silva et al. 1998) normalized to the 850/xm point and non-Spitzer data from the literature (upper limits at 3-cr). 
Both sources have broad-band SEDs which are similar to Arp 220. Dotted lines show a simple power law, f u ocv a , fit to the 3.6-8/xm data. We 
also show the best fit (stellar) SED from HYPERZ, fitted to the optical/mid-infrared photometry (see §4 for details). 



z = 3.1, it is possible that LAB la and b are interacting given 
their projected separation of just ~4" (~30kpc at z = 3.1). 

Previous multi-wavelength observations of LAB 1 (Chap- 
man et al. 2004) associated the SMG with the brightest of 
a group of UV fragments near the geometrical center of the 
structure and close to the 21 cm source position (Fig. 1). Com- 
ponent LAB la is in the vicinity of this group of faint UV 
fragments, but the 3.6-8/im emission is coincident with a 
Ks = 21.5 Extremely Red Object (ERO) identified in Steidel 
et al. (2000), rather than the brightest UV component ('Jl' 
in Chapman et al. 2004) which is closest to the radio peak 
(and a tentative CO(4-3) detection). Both LAB la and the 
UV object Jl are plausible sources of the submm emission. 
However, given that SMGs are expected to be massive, we 
view the 3.6-8/im detection as a more reliable indicator of 
the SMG's counterpart than the weak UV source. We there- 
fore propose that LAB la is the most likely IRAC counterpart 
to the SMG in LAB 1 

There are three mid-infrared counterparts within the extent 
of LAB 2 (Fig. 1, Table 1). We can eliminate LAB 2c imme- 
diately as a potential counterpart given its blue mid-infrared 
colors and large, bright optical counterpart in the HST image 
- it is a foreground object. The remaining two counterparts 
LAB 2a and LAB 2b have similar 3.6-5.8/im characteristics, 
and similar compact HST morphologies (Fig. 1). However, 
only LAB 2b is detected at 8/im, which has been shown to be 
a useful identifier of the mid-infrared counterparts of SMGs 
(Ashby et al. 2006). Moreover, LAB 2 contains a hard X-ray 
source in the halo (Basu-Zych & Scharf 2004), with an esti- 
mated luminosity L x ~ 10 44 erg s -1 . The location of the X-ray 
source is consistent with LAB 2b (AO = 0.8 ± 1.4"), whereas 
LAB 2a is in the extreme south-east of the LAB, ~8.5" away 
from the X-ray source. Given the association between SMGs 
and moderately bright X-ray sources (Alexander et al. 2005a), 
the SMG in LAB 2 is likely to lie close to the position of the 
X-ray source. We therefore propose that LAB 2b is the mid- 
infrared counterpart to the SMG in LAB 2 8 . The precise posi- 

8 LAB 2b is also at the location of the LBG M 14, although in this case 
the LBG appears kinematically associated with the wider halo (Wilman et al. 



tions of the SMGs in LAB 1 and LAB 2 and the offsets to the 
other mid-infrared components are given in Table 1 . 

4. CLASSIFICATION OF THE SMGS 

We can use the 3.6-8/im photometry (along with multi- 
wavelength data for these objects from the literature; Chap- 
man et al. 2001, 2004; Steidel et al. 2000) to construct the 
optical-radio SEDs for the SMGs (Figure 2). We present 
new / and K band photometry from the UKIDSS-Deep Ex- 
tragalactic Survey (to be discussed in Stott et al. 2007 in 
prep). We note that the poor constraints on the /24//850 ratio 
means that we cannot make robust conclusions about the rest- 
frame mid-to-far infrared differences between LAB la and 
LAB 2b. Both have /24//850 upper limits consistent with ob- 
scured ULIRGs at z ~ 3 (Pope et al. 2006). Similarly the 
limits of the observed fs.s/fsso ratios for LAB la and LAB 2b 
are consistent with other SMGs at z ~ 3, although the ratio for 
LAB la is at the low end of the distribution (Pope et al. 2006) 
indicating either a relatively low stellar mass for this system, 
or more likely very high obscuration. Thus, the SMGs asso- 
ciated with LABs are similar to the wider SMG population. 

We now examine the rest- frame near-infrared portion of the 
SEDs using the IRAC data. Near-infrared colors have been 
shown to be a useful diagnostic for the power-source of local 
luminous infrared galaxies (e.g. Spinoglio et al. 1995). For 
LAB la we find rest- frame (Vega) colors (J-H) = 0.6 ± 0.5 
and (H-K) = 0.8 ±0.4, and in LAB 2b we find (J-H) = 0.8 ± 
0.6 and (H-K) = 1 . 1 ± 0.4. LAB la's colors are similar to that 
of 'normal' galaxies ((J-H) ~ 0.8, (H-K) ~ 0.3; Spinoglio 
et al. 1995), whereas LAB 2b exhibits a redder rest- frame (H- 
K) color, indicating an excess in the rest-frame ^-band likely 
associated with a non-stellar contribution. We fit a simple 
power-law, f v ex v a , to the IRAC photometry for LAB 2b, 
obtaining a = -1.0±0.2, similar to that expected for AGN 
(Alonso-Herrero et al. 2006). 

2005). We note that in the submillimeter observation of this halo, LAB 2b is 
actually situated near to the edge of the SCUBA beam (Fig. 1). The correction 
to the 850^m flux to account for this offset would make the SMG brighter, 
with ~ 5 mJy (Chapman et al. 2001). 
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MlD-IR PROPERTIES OF SOURCES LYING WITHIN LABS 
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NOTE. — (1) Positional offsets from the SMGs in LAB 1 
and LAB 2 identified in §3. The proposed mid-infrared counter- 
parts, LAB la and LAB 2b are at 22 h 17 m 26.0 s , +00°12 / 36.2 // and 
22 h 17 m 39.1 s , +00°13 / 30.4 // (J 2000) respectively. Positional uncer- 
tainties based on the 3.6/xm centroids are 0.8". (2) Flux uncertain- 
ties are 0.3/xJy in 3.6/xm & 4.5/xm and 2.2/iJy in 5.8/xm & 8.0/xm. 
Quoted limits are 3 -a. 

To further test if the optical-mid-infrared SEDs of the 
sources can be reproduced by stellar populations, we use 
hyperz (Bolzonella, Miralles & Pello 2000) to model the 
stellar continuum by fitting to the rest- frame 0.15-2/im pho- 
tometry. We adopt a model with an exponentially declin- 
ing SFR, with r = 30Gyr, from Bruzual & Chariot (1993) 
and constrain the redshift to be z = 3.09. Dust is modelled 
with a Calzetti law (Calzetti et al. 2000). The rest- frame 
optical-near-infrared portions of the SEDs are well fit by this 
model, with high-reddening (A v ~ 3) and young ages. Ab- 
solute K-bznd magnitudes based on the fit are M K ~ -26.8 
in LAB la and M K ~ -27.2 for LAB 2b. These absolute rest- 
frame AT-band magnitudes are comparable to the wider SMG 
population (Borys et al. 2005), and suggest stellar masses of 
~ 10 11 M , although the /T-band magnitude of LAB 2b is 
likely to be contaminated by AGN emission, with a significant 
rest- frame 2/im excess over the hyperz fit of ~5/iJy. The 



detection of X-ray emission from LAB 2b is the best evidence 
we have that it contains an AGN. To quantitatively test the 
correspondence between the mid-infrared excess and the X- 
ray emission we use the rest-frame X-ray-10.5/im correlation 
for AGN from Krabbe et al. (2002). We predict the rest-frame 
10.5/im emission for LAB 2b by extrapolating the rest- frame 
2//m using the measured spectral slope of a = -1. With this 
approach we estimate the X-ray luminosity from the AGN to 
be of order 10 44 erg s -1 , in reasonable agreement with that ob- 
served (Basu-Zych & Scharf 2004). The luminosity of the 
X-ray source is consistent with a moderate contribution to the 
bolometric emission from an AGN (Alexander et al. 2005b), 
however it is still not clear if an AGN or starburst is responsi- 
ble for powering the extended Lyman-a emission. LAB la's 
rest-frame AT-band excess over the best fit stellar template, of 
~3±2/iJy is less significant. Nevertheless, taking the same 
approach as for LAB 2b, we predict an X-ray luminosity of 
~5 x 10 43 ergs _1 , consistent with the 2-10 keV limit in the 
existing Chandra observations. 

To make progress we need mid-infrared observations of a 
larger sample of LABs in SA22, four of which also contain 
SMGs (Geach et al. 2005). This will allow us to investigate 
whether they also exhibit rest- frame near-infrared characteris- 
tics similar to the general submillimeter population. Searches 
for obscured AGN in SA22 will also be significantly ad- 
vanced by a recently approved 400 ks Chandra ACIS-I ex- 
posure of this region. Together these studies will provide a 
clearer view of the energetics of LABs and their relationship 
to the wider SMG population. 
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